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ABSTRACT 

The  Department  of  Defense  is  concerned  with  the  viability 
of  an  airborne  high  energy  laser  system.   The  laser  is  housed 
in  a  blunt  turret  atop  a  NKC-135  aircraft.   Turbulence  gen- 
erated by  flow  separation  around  the  turret  causes  optical 
distortion  of  the  laser  beam.   Control  of  flow  separation  is 
needed  to  improve  laser  beam  performance  especially  for  aft- 
aimed  turrets . 

One  technique  proposed  for  flow  control  is  a  fairing  de- 
sign which  will  stabilize  shed  vortices  by  suction.   A  two 
dimensional  computer  model  was  used  to  design  a  fairing  com- 
patible with  present  test  equipment. 

Experimental  research  of  this  fairing  design  was  con- 
ducted in  wind  tunnel  tests.   Although  flow  mapping  demon- 
strated improved  flow  performance  through  the  use  of  suction, 
total  quiescent  flow  was  never  achieved.   A  more  adequate 
three  dimensional  model  is  needed  to  design  a  fairing  that 
will  stabilize  trapped  vortices. 


COMMENT  CONCERNING  JOINT  RESEARCH  EFFORT 

This  thesis  and  LCDR  David  Rippel's  thesis,  Airborne 
Laser  Turret  Flow  Control:   A  Parametric  Study  of  Wind  Tun- 
nel Conditions  [Ref.  1]  ,  were  the  result  of  a  joint  research 
effort.   The  basic  data  acquisition  system,  was  the  same  for 
both  projects,  however  this  thesis  deals  with  the  quality  of 
flow  using  a  nosepiece  designed  to  stabilize  shed  vortices. 
The  results  in  Ref.  1  are  a  more  complete  testing  of  two 
earlier  nosepiece  designs. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Over  the  years  the  Department  of  Defense  (DOD)  has  been 
concerned  with  the  deployment  of  a  high  energy  laser  as  a 
viable  weapons  system.   Airborne  laser  research  is  being 
conducted  at  the  Air  Force  Weapons  Laboratory,  Kirtland  Air 
Force  Base,  New  Mexico.   Two  technically  complex  NKC-135 
aircraft  are  being  used  as  the  airborne  laboratories. 

The  high  energy  laser  weapon  system  radiates  large 
amounts  of  thermal  energy  on  another  airborne  target  ren- 
dering it  useless.   The  main  components  of  the  system  include 
the  laser  itself  and  the  beam  control  subsystem.   The  turret 
on  top  of  the  NKC-135  aircraft  houses  the  pointer-tracker 
assembly  which  aims  the  laser  beam.   Due  to  the  blunt  body 
design  inherent  with  the  rounded  laser  turret,  shear  layers, 
boundary  layers,  flow  separation,  and  vortex  shedding  can 
cause  noticeable  distortions  in  the  laser  beam  energy. 
Unsteady  pressure  loads  on  the  turret  and  optical  components 
cause  jitter,  a  phenomenon  that  translates  beam  energy  fo- 
cused on  a  small  area  into  a  larger  area.   The  time  needed 
to  inflict  damage  to  the  target  is  increased. 

Separated  flow  must  be  controlled  in  order  to  enhance 
optical  qualities  of  the  laser  beam.   Research  and  experi- 
mentation have  demonstrated  that  optical  distortion  caused 
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by  unsteady  flow  cannot  be  corrected  by  adaptive  optical 
systems  due  to  the  bandwidth  requirements  exceeding  current 
technology. 

Various  methods  for  attempting  flow  control  over  a  tur- 
ret were  introduced  at  a  workshop  sponsored  by  Captain 
Richard  deJonckheere  at  the  Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  New  Mexico  [Ref.  2].   One  innovative 
method  discussed  in  this  workshop  was  an  attempt  to  trap 
vortices  shed  by  the  turret  and  stabilize  them  by  suction 
through  ports  located  in  the  most  optimum  locations.   This 
stabilization  could  then  hopefully  create  quiescent  flow 
patterns  over  the  surface  of  the  turret.   Figure  1-1  is  the 
original  design  submitted  by  Craig  [Ref.  3], 

B.   THESIS  OBJECTIVE 

The  primary  objective  of  this  thesis  research  was  to 
design,  manufacture,  and  test  a  f airing-nosepiece  device 
which  would  employ  suction  to  stabilize  trapped  vortices 
shed  by  a  laser  turret.   A  design  compatible  with  previous 
test  equipment  was  necessary.   The  suction  holes  would  be  in 
the  wind  tunnel  floor. 
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II.   VORTEX  THEORY 

A.   VORTEX  GENERATION 

1 .   Boundary  Layer  Separation 

Viscosity  is  that  property  of  a  real  fluid  which 
creates  shear  forces  between  two  fluid  elements.   Air  and 
water  both  have  small  viscosities;  therefore  in  earlier  days, 
scientists  assumed  that  forces  due  to  viscous  friction  were 
negligible  compared  with  pressure  and  gravity  forces.   The 
science  of  classical  hydrodynamics  is  based  on  the  premise 
that  fluid  is  ideal,  i.e.,  frictionless  and  incompressible. 
Empirical  results,  however,  did  not  always  agree  with  clas- 
sical fluid  theories.   This  is  especially  true  when  a  solid 
body  is  immersed  in  a  fluid  flow.   Classical  hydrodynamic 
calculations  lead  to  zero  drag  on  the  immersed  body,  which 
is  impossible.   This  contradiction  is  referred  to  as  d'Alem- 
bert ' s  paradox,  and  it  was  two  centuries  before  a  solution 
to  the  problem  was  discovered.   The  concept  of  the  boundary 
layer,  theorized  by  L.  Prandtl  in  1904  [Ref .  4] ,  bridged 
this  scientific  gap.   He  deduced  that  the  fluid  flow  is  di- 
vided into  two  regions.   The  fluid  near  the  solid  surface  is 
in  a  thin  boundary  layer  where  the  velocity  gradient  is 
large  enough  that  friction  within  the  boundary  layer  due  to 
fluid  viscosity  must  be  considered.   The  other  region  is 
outside  the  boundary  layer  and  the  effect  of  fluid  viscosity 
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is  negligible.   Fluid  in  the  outer  region  follows  classical 
potential  flow  theories. 

Pressure  is  a  force  exerted  on  a  unit  area  and  a 
pressure  gradient  is  a  change  in  pressure  along  a  particular 
direction.   When  there  is  an  adverse  pressure  gradient  (pres- 
sure increasing  in  the  downstream  direction)  the  boundary 
layer  thickens  rapidly  along  the  surface.   Adverse  pressure 
gradients  act  to  oppose  the  downstream  fluid  motion  and 
reduce  the  fluid  momentum  in  the  boundary  layer  area.   High 
viscous  frictional  forces  in  the  area  next  to  the  surface 
tend  to  retard  the  flow  velocity  even  more.   The  combination 
of  the  opposing  pressure  and  the  viscous  frictional  forces 
can  ultimately  cause  the  phenomenon  of  flow  separation. 
When  the  flow  separates,  large  turbulent  eddies  can  be 
formed . 

The  position  of  the  point  of  separation  is  dependent 
on  the  roughness  of  the  body,  the  Reynolds  number  of  the 
flow,  and  the  geometric  curvature  of  the  body.   By  examining 
an  irrotational  form  of  Bernoulli's  equation 


Pl    Vl2    p2    V22 

—  +  -£-  =  —  +  -4-  (2-1) 

P1  2  P2     2 


it  can  be  seen  that  adverse  pressure  gradients  are  generally 
created  when  there  is  a  deceleration  of  the  fluid  motion. 
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For  example,  if  V2  is  smaller  than  V,,  then  p2  is  larger 
than  p1 ,  providing  the  density  does  not  change. 

Flow  mapping  over  a  cylinder  can  be  seen  in  Figure 
II-l,  taken  from  Pao  [Ref.  5].   Part  (a)  shows  the  stream- 
line pattern  and  pressure  distribution  for  a  non-viscous 
flow.   Part  (b)  shows  a  laminar  boundary  layer  separation 
due  to  the  decelerated  flow  in  the  rear  portion  of  the  cyl- 
inder.  Part  (c)  also  shows  a  separated  flow,  but  the  sepa- 
ration points  are  farther  downstream  caused  by  an  increase 
in  fluid  momentum  generated  by  a  turbulent  flow. 
2 .   Von  Karman  Vortex  Street 

Von  Karman  noted  that  when  a  cylindrical  body  is 
moved  through  a  fluid,  a  series  of  alternating  vortices  are 
generated.   Generally  the  vortices,  along  one  line  all  have 
circulation  in  the  same  direction,  but  are  opposite  to  the 
vortices  in  the  other  line  [Ref.  6].   A.  Roshko  at  California 
Institute  of  Technology,  Pasadena,  also  studied  the  vortex 
street  phenomenon  and  determined  that  the  formation  and 
shape  of  the  vortices  are  highly  dependent  on  the  Reynolds 
number  [Ref.  7,  8].   At  higher  Reynolds  numbers,  turbulence 
breaks  down  vortex  street  periodicity.   During  this  process 
the  distance  between  individual  vortices  increases  steadily 
while  the  circulation  decreases  [Ref.  8:   p.  222]. 
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B.   VORTEX  MOTION 

1.  Forced  Vortex  Motion 

A  forced  vortex  motion  occurs  when  a  constant  torque 
is  maintained  on  a  body  of  fluid.   This  is  illustrated  when 
a  circular  container  of  fluid  is  rotated  about  its  central 
axis  with  a  constant  angular  momentum  [Ref.  5:   pp.  139-140] 
The  velocity  of  any  fluid  particle  is  proportional  to  the 
radius  outward  from  the  axis  of  rotation 

v  =  rC,  (2-2) 

Figure  II-2  shows  the  velocity  distribution  of  a  forced 
vortex  motion. 

2 .  Free  Vortex  Motion 

A  free  vortex  motion  occurs  naturally  when  fluid 
masses  move  in  curved  paths.   The  velocity  of  a  free  vortex 
motion  is  seen  to  vary  inversely  with  the  radial  distance 
from  the  vortex  center 


C2 
V  =  ~  (2-3) 


Figure  II-3  depicts  the  free  vortex  velocity  distribution. 
3 .   Rankine  Combined  Vortex 

The  velocity  of  flow  in  a  free  vortex  would  go  to 
infinity  as  the  radius  goes  to  zero.   This,  however,  is  phy- 
sically impossible  due  to  viscous  forces  at  the  central  re- 
gion.  The  resulting  flow  resembles  the  simple  case  Rankine 
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combined  vortex,  which  is  a  combination  of  forced  and  free 
vortex  motion.   Circular  air  flow  in  a  tornado  is  an  ap- 
proximation of  a  Rankine  combined  vortex.   Figure  II-4  is 
an  illustration  of  the  Rankine  combined  vortex  velocity 
distribution. 

C.   IMPROVED  VORTEX  VELOCITY  DISTRIBUTION 
1 .   Core  Flow 

Various  efforts  have  been  described  in  Lewellen 
[Ref.  9]  to  improve  upon  the  velocity  distribution  in  order 
to  smooth  the  sharp  discontinuity  between  the  combined 
forced  and  free  vortex  motions.   The  Rankine  vortex  is  the 
simplest  solution  to  the  angular  momentum  equation  noted  in 
Ref.  9.   The  angular  momentum  equation  is 


P(vr>      fa   1  3(vr)    £>l 

Dt     Jr  [3r   r   3r   +  ~7 

1  3z 


This  equation  does  assume  that  the  core  of  the  tornado  is 
axisymmetric  and  that  the  eddy  viscosity,  v,  is  constant. 
2 .   Burgers-Rott  Solution 

Lewellen  [Ref.  9]  notes  that  a  well-behaved  solution 
to  equation  (2-4)  was  obtained  by  Burgers  and"  Rott  independ- 
ently.  For  their  solution,  consider  a  cylindrical  coordinate 
system  with  velocity  components  u,  v,  and  w.   The  radial  in- 
flow velocity,  u,  is  defined  as 

u  =  -ar  (2-5) 
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The  axial  flow  velocity,  w,  is 

w  =  2az  (2-6) 

while  the  tangential  velocity,  v,  is  defined  as 


r 

00 

r 


i  -ar 
1  -  exp 


2v  /J 


(2-7) 


The  variable,  a,  used  in  all  three  velocity  components,  is 
the  upflow  gradient  and  has  the  units  of  1/sec.  The  pres- 
sure distribution  consistent  with  this  solution  is 


2 
ar 

p(r,z)  =  p(0,0)  +  p-^— J     (1-^   )  dx 

0 


2 
£§-(r2  +  4z2)  -  pgz  (2-8) 


The  kinematic  viscosity,  v,  in  equation  (2-7)  becomes  the 
effective  turbulent  kinematic  viscosity  when  the  flow  in  the 
vortex  is  turbulent.   A  characteristic  velocity  for  the  flow 


is  /a/v  r  and  a  characteristic  length  is  vv/a.  The  velocity 
and  radial  distance  in  Figure  II-5  have  been  normalized  with 
respect  to  these  characteristic  quantities. 
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III.   NOSEPIECE  DESIGN 

A.   DESIGN  PARAMETERS 

Refer  to  Figures  III-l  and  III-2,  where  the  four  compo- 
nents are  identified.   These  components  include:   turret, 
nosepiece,  fairing,  and  plenum  for  suction  at  wind  tunnel 
floor. 

After  the  initial  decision  was  made  to  study  the  possi- 
bilities of  the  trapped  vortex  nosepiece,  a  design  process 
was  begun.   The  new  nosepiece  was  to  be  completely  compati- 
ble with  turret,  fairing,  and  blower  dimensions  in  order  to 
prevent  major  modifications  to  the  existing  system.   Results 
of  the  trapped  vortex  research  could  also  be  compared  more 
easily  with  data  obtained  during  experimentation  of  the  Ta- 
pered Symmetric  Nosepiece  (TSN)  and  the  Uniform  Conformal 
Nosepiece  (UCN)  [Ref.  1,  11,  and  12]. 

Fixed  apparatus  measurements  included  a  fairing  width  of 
10.5  inches,  a  turret  diameter  of  16.8  inches,  and  a  plenum 
area  18  by  20  inches.   A  research  effort  was  dedicated  to 
determine  critical  values  for  the  geometry.   Figure  III-l  is 
a  top  view,  looking  down  on  the  wind  tunnel  floor.   The  solid 
lines  represent  fixed  apparatus  measurements,  while  the 
dashed  lines  represent  variable  geometric  parameters.   The 
length,  L,  is  the  distance  between  the  center  of  the  turret 
and  the  fairing.   The  fairing  was  fixed  to  the  tunnel  floor 
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and  could  not  be  moved.   Nosepiece  dimensions  are  determined 
by  the  length,  L- ,  and  the  width,  W.   Suction  port  locations 
are  determined  by  R~  and  W- .   The  variable,  D,  is  the  suc- 
tion port  diameter,  while  R  is  the  radius  that  would  define 
the  curvature  of  the  nosepiece.   The  turret-nosepiece- 
f airing  assembly  can  also  be  seen  in  Figure  III-2.   Note 
that  the  values  L^  and  L.    determine  the  nosepiece  length. 
Top  nosepiece  curvature  is  represented  by  the  value  R., , 
which  is  at  a  height,  H,  above  the  tunnel  floor.   Suction 
port  area,  vortex  location,  and  nosepiece  curvature  were 
three  of  the  more  important  design  considerations. 

B.   ELECTROSTATIC  DESIGN 

During  the  early  design  stages,  a  method  was  needed  to 
visualize  flow  patterns  as  influenced  by  different  geometric 
design  parameters.   Electrostatic  plotting  seemed  like  an 
effective  tool  to  use  for  problem  modeling.   The  PASCO  Model 
224  Electrometer  manufactured  by  PASCO  Scientific,  San  Lean- 
dro ,  California,  was  used  as  a  flow  plotting  device.   Various 
nosepiece  and  fairing  designs  were  tested  in  order  to  deter- 
mine the  effectiveness  of  this  system.   Figure  III-3  shows 
the  arrangement  of  the  test  design  between  the  two  elec- 
trodes.  The  electrodes  and  the  test  design  were  drawn  with 
a  silver  conducting  pen.   Two  leads  were  joined  to  points  A 
and  B  in  Figure  III-3  and  were  connected  to  a  variable  poten- 
tiometer.  Increasing  current  flow  in  these  leads  would  be 
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akin  to  increased  suction  through  base  platform  vortex 
ports.   Lines  of  equipotential  correspond  to  streamlines. 
The  lines  of  equipotential  which  were  then  plotted  around 
the  entire  test  design  represent  the  streamlines  around  an 
aerodynamic  surface.   However,  the  test  apparatus  was  not 
precise  enough  in  the  area  of  the  nosepiece  to  make  design 
decisions.   Further  research  attempts  using  the  electro- 
static analogy  were  cancelled. 

C.   COMPUTER  DESIGN 

Lack  of  success  in  the  electrostatic  design  phase  and 
the  complexity  of  three  dimensional  mathematical  models 
forced  the  use  of  a  computer  program  geared  to  two  dimen- 
sional flow. 

1.   Panel  Method 

A  Fortran  program,  which  was  developed  by  Adjunct 
Professor  Peter  Bellamy-Knights,  applied  the  panel  method 
to  determine  velocity  distributions  and  pressure  coeffi- 
cients over  any  aerodynamic  surface.   The  surface  could  be 
divided  into  panels,  each  with  its  own  source  strength.   An 
iterative  routine  using  linear  equations  in  matrix  format 
solved  for  the  magnitude  of  the  source  strengths,  the  velo- 
cities tangential  to  the  panels,  and  the  corresponding  pres- 
sure coefficients. 

In  order  to  use  the  panel  method  for  design,  various 
additions  were  incorporated.   The  turret-nosepiece-f airing 
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assembly  was  assumed  to  be  a  one-piece  body,  symmetrical 
with  respect  to  the  X-axis.   Two  vortices  of  variable  posi- 
tion and  strength  were  located  outside  the  nosepiece  area; 
vortex  location  is  shown  in  the  top  half  of  Figure  III-4. 
The  effects  of  these  two  vortices  influenced  the  magnitude 
of  all  the  other  panel  strengths,  which  in  turn  determined 
the  final  output.   Design  variations  were  necessary  in  the 
region  of  the  nosepiece,  therefore  different  coordinates 
could  be  input  to  describe  the  different  panels.   By  varying 
nosepiece  geometry,  it  was  possible  to  deduce  flow  trends 
along  the  surface  of  the  panels. 

The  two-dimensional  flow  conditions  calculated  using 
the  panel  method  are  not  a  precise  representation  of  the 
actual  flow  for  two  reasons.   First,  the  actual  geometry  is 
three-dimensional,  and  second,  the  calculation  ignores  vis- 
cous effects.   Recognizing  these  deficiencies,  a  design 
philosophy  was  used.   The  design  philosophy  is  illustrated 
in  Figure  III-4.   Refer  to  the  plane  of  symmetry  in  Figure 
III-4.   Below  the  plane  of  symmetry  is  an  equivalent  solid 
body  which  results  from  a  trapped  vortex.   Above  the  plane 
of  symmetry  is  the  flow  with  the  trapped  vortex.   The  line 
segment  A-B  in  Figure  III-4  is  the  dividing  streamline  be- 
tween the  external  flow  and  the  trapped  vortex.   Along 
streamline  A-B  is  a  distribution  of  pressure  and  velocity. 
The  aim  of  the  computer  model  was  to  match  flow  conditions, 
i.e.  velocity  and  pressure,  along  streamline  A-B.   The  flow 
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conditions  along  A-B  due  to  the  trapped  vortex  should  be 
identical  to  flow  conditions  along  A'-B*  due  to  an  equiva- 
lent solid  body.   A  copy  of  the  Fortran  computer  program  is 
included  at  the  end  of  this  report. 
2 .   Input  Parameters 

In  order  to  allow  easy  transfer  of  coordinate  inputs, 
all  values  were  non-dimensionalized  with  respect  to  the  tur- 
ret radius.   The  turret-nosepiece-f airing  assembly  was  di- 
vided into  thirty  total  panels.   More  panels  were  positioned 
in  areas  of  rapid  slope  changes  such  as  the  nosepiece  and 
turret  curves.   Figure  III-5  shows  the  basic  breakdown  of 
the  initial  panel  design.   Table  III-l  gives  the  x  and  y 
coordinate  values  of  all  thirty  points.   Points  6  through  12 
and  their  negative  counterparts,  20  through  26,  are  com- 
pletely variable  and  determine  a  design  geometry  of  the 
nosepiece.   Figures  III-6,  III-7  and  III-8  illustrate  other 
nosepiece  variations  which  were  calculated.   Table  III-2 
gives  the  coordinate  locations  of  these  three  other  designs. 

A  model  was  needed  to  describe  vortex  strength. 
Consider  fluid  flow  along  a  line  in  a  volume.   The  dot  pro- 
duct of  the  instantaneous  velocity  vector  and  a  line  element 
will  produce  the  component  of  the  velocity  along  that  line. 
Circulation  is  the  line  integral  of  that  dot  product  along  a 
closed  path.   The  mathematical  expression  for  circulation  is 


r  =jfu«ds  (3-1) 
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where  U  is  the  instantaneous  velocity  and  ds  is  the  line 
element.   This  equation  can  be  non-dimensionalized  by  divi- 
ding both  sides  by  the  freestream  velocity,  Uao,  and  a 
characteristic  length.   The  characteristic  length  in  this 
case  is  the  turret  radius,  R.   The  non-dimensional  equation 
becomes 


*  =Md(t 


(3-2) 


For  a  free  vortex,  the  circulation  is  proportional 
to  the  strength  of  the  free  vortex.   The  computer  program 
uses  a  variable  F,  which  is  equal  to  the  value  T/U^R  from 
equation  (3-2).   Thus,  F  is  used  to  describe  non-dimensional 
vortex  strength.   As  noted  earlier,  the  two  vortex  locations 
and  strengths  were  completely  variable.   The  computer  pro- 
gram inputs  vortex  data  as  X2 ,  Y2 ,  and  F  as  the  X-coordinate , 
Y-coordinate,  and  non-dimensional  vortex  strength  respec- 
tively.  The  negative  counterpart  is  automatically  calculated 
in  the  computer  program.   The  value  of  (F)  was  assumed  to  be 
0.6  and  kept  as  a  constant  parameter. 
3 .   Output  Parameters 

The  most  important  output  parameter  from  the  program 
was  the  tangential  velocity  vectors  produced  along  the  wall 
of  the  nosepiece  panels.   Plots  were  made  of  these  vectors 
to  determine  the  direction  and  strength  of  the  flow  field. 
The  freestream  velocity  was  taken  as  unity,  thus  a  velocity 
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downstream  could  be  compared  to  it.   For  example,  if  a  velo- 
city vector  read  a  value  of  -0.5,  the  tangential  velocity  of 
that  panel  at  its  midpoint  would  be  one-half  the  freestream 
value  and  in  the  opposite  direction. 

Figure  III-7  demonstrated  better  flow  conditions 
during  calculations,  but  the  design  needed  to  be  shortened 
to  produce  a  stronger  region  of  reversed  flow  along  the 
panels . 

Figure  III-9  shows  a  plot  of  the  final  configuration 
and  the  tangential  velocity  vectors  determined  from  the  com- 
puter output.   This  configuration,  coupled  with  the  vortex 
location  specified,  proved  to  be  the  only  design  with  a 
velocity  magnitude  of  zero  along  the  panel  between  points 
eleven  and  twelve.   All  other  designs  tested  had  a  larger 
flow  vector  in  the  direction  towards  point  twelve.   Note 
how  the  velocity  vectors  changed  direction  at  point  A  in 
Figure  III-9.   It  was  assumed  that  the  flow  would  be  forced 
to  travel  along  a  path,  similar  to  the  line  A-B.   This  flow 
pattern  is  typical  of  the  design  philosophy  noted  before  in 
Figure  III-4. 

The  coordinates  of  the  selected  design  were  faired 
in,  and  a  three  dimensional  body  was  constructed.   The  for- 
ward portion  of  the  nosepiece  was  cut  to  meet  the  aft  por- 
tion of  the  turret.   Figures  111-10,  III-ll,  and  111-12  show 
the  final  completed  side,  top,  and  end  views  respectively. 
Table  III-3  gives  the  coordinates  of  the  30  points  used  in 
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the  final  panel  design,  while  Table  III-4  gives  the  tangen- 
tial velocities  and  pressure  coefficients  along  the  panels 
of  the  design  shown  in  Figure  III-9. 
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IV.   EXPERIMENTAL  APPARATUS 

A.  WIND  TUNNEL 

Wind  tunnel  tests  were  conducted  in  the  Naval  Postgradu- 
ate School  five-foot  by  five-foot  low-speed  tunnel  at  a 
maximum  velocity  of  33  feet  per  second.   Using  the  pre- 
viously built  one-third  scale  turret  model  (D  =  16.8  inches) 
a  Reynolds  number  of  about  3  x  10   was  achieved.   According 
to  Schlichting  [Ref.  10:   pp.  20-21],  the  value  of  the  Rey- 
nolds number  for  the  tests  was  in  the  critical  range  and 
flow  should  be  turbulent  in  the  boundary  layer  on  the  turret. 

B .  BLOWER 

The  blower  was  chosen  during  previous  thesis  research 
[Ref.  11  and  12],  and  is  the  Backward  Inclined  Airfoil,  Model 
500,  Single  Width  Single  Inlet  (B.I. A. -500,  SWSI)  made  by  the 
Aerovent  Company  Inc.,  of  Piqua,  Ohio.   This  centrifugal 
blower  has  a  capacity  of  7700  cfm  with  a  static  pressure 
differential  of  14  inches  H-0.   A  control  valve  on  the  blower 
face  could  be  operated  from  a  full  closed  to  full  open  posi- 
tion, allowing  for  an  increase  in  airflow.   A  sheet  metal 
structure  was  used  to  mate  the  blower  inlet  to  the  wooden 
ductwork  underneath  the  tunnel  test  section.   Figure  IV- 1 
shows  the  front  side  of  the  blower  and  completed  ductwork 
below  the  wind  tunnel. 
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C.  DUCTING 

A  wooden  rectangular  duct  connected  the  hollow  fairing 
to  the  blower  assembly.   Four  adjustable  butterfly  valves 
were  positioned  inside  four  separate  flow  channels  within 
the  ductwork.   A  fifth  channel  led  to  a  large  plenum  chamber 
at  the  base  of  the  turret  which  was  used  in  previous  research 
as  a  path  for  fuselage  boundary  layer  bleed.   Figure  IV- 2 
shows  the  fairing  positioned  behind  the  plenum  chamber  with 
the  nosepiece  and  turret  removed.   The  plenum  duct  also  had 
a  fifth  butterfly  valve  which  was  locked  in  the  full  open 
position.   Figure  IV- 3  shows  the  back  side  of  the  ductwork 
under  the  tunnel  flooring  and  the  five  adjustable  butterfly 
valves . 

D.  FAIRING 

A  filler  piece  six  inches  long  was  used  to  fill  the  gap 
between  the  previous  fairing  position  and  plenum  cavity  on 
the  tunnel  floor.   The  plenum  chamber  could  now  be  used  as  a 
vortex  suction  unit  without  major  modification  of  existing 
ductwork.   Figure  IV- 4  shows  the  filler  piece  construction. 
The  hollow  fairing  was  not  used  as  a  suction  device,  so  the 
four  other  butterfly  valves  were  locked  in  the  fully-closed 
position.   The  fairing  was  10.5  inches  wide,  which  produced 
a  maximum  look-back  angle  of  150  degrees. 
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E.  TURRET 

A  hollow,  one-third  scale  model  of  the  existing  airborne 
laser  turret  was  constructed  on  drawings  provided  by  Captain 
Richard  de Jonckheere .   The  wooden  model  consists  of  a  hollow 
16.8  inch  diameter  circular  cylinder,  9.6  inches  in  height, 
topped  by  a  16.8  inch  diameter  hemisphere.   Twenty-five  in- 
dividual pressure  taps  were  placed  over  the  area  of  the 
turret. 

F.  NOSEPIECE 

The  nosepiece  was  constructed  from  sheet  metal,  wood,  and 
styrofoam.   It  was  placed  in  between  the  fairing  filler  sec- 
tion and  the  laser  turret.   The  upper  one-third  of  the  model 
was  a  three  dimensional  body  based  on  the  design  criteria 
from  Chapter  III.   Another  nine  pressure  taps  were  located 
on  the  left  side  of  the  model.   Figures  IV-5  and  IV-6  show 
the  general  shape  of  the  nosepiece. 

G.  INSTALLATION 

The  plenum  chamber  was  covered  by  a  27  x  27  inch  square 
sheet  of  one-half  inch  plywood  in  order  to  create  a  base 
platform  for  the  nosepiece  and  turret.   Figures  IV-7  and 
IV-8  show  various  views  of  the  fairing,  filler  piece,  turret 
and  nosepiece.   Figure  IV-9  shows  a  side  view  drawing  of  the 
entire  test  setup.   Various  size  holes  were  drilled  in  the 
base  platform  in  order  to  suck  the  trapped  vortices  using 
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the  blower  assembly.   Figure  IV- 10  shows  the  final  configura- 
tion of  the  base  platform  and  the  adjustable  hole  control. 
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V.   INSTRUMENTATION 

A.  PRESSURE  TAPS 

Pressure  taps  were  installed  in  the  turret,  wind  tunnel, 
and  nosepiece.   Due  to  the  varied  location  of  all  the  pres- 
sure taps,  a  pressure  distribution  of  the  turret  and  nose- 
piece  surfaces  could  be  plotted.   Figure  V-l  is  a  pictorial 
representation  of  the  various  pressure  tap  locations  on  the 
turret  surface,  while  Figure  V-2  shows  the  pressure  tap  lo- 
cations on  the  left  side  of  the  nosepiece.   All  pertinent 
pressure  taps  were  connected  directly  to  a  48  port  Scani- 
valve  via  flexible  Tygon  tubing. 

B.  DATA  ACQUISITION  SYSTEM 

The  wind  tunnel  data  acquisition  system  used  during  the 
course  of  research  included  an  INTEL  80/10  Computer  System, 
an  AN/UGC-59A  Teletypewriter  Set,  a  48  port  Scanivalve ,  and 
an  analog  to  digital  display  box.   Figure  V-3  shows  the  com- 
puter system,  teletypewriter,  and  digital  display  unit. 
Figure  V-4  is  a  picture  of  the  Scanivalve. 

A  master  control  program  for  Scanivalve  operation  was 
developed  to  ensure  that  voltage  readings  taken  from  the 
system  could  be  displayed  on  the  teletypewriter.   These 
voltage  readings  from  a  capacitor  pressure  transducer  in 
the  Scanivalve  were  converted  to  pressure  in  centimeters  of 
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water  using  a  U-tube  manometer.   Due  to  the  linearity  of 
the  various  voltages  plotted  against  different  pressure 
values  of  the  U-tube,  an  equation  was  determined  to  facili- 
tate conversion 

y  =  9.2608x  +  0.0269 

where  x  is  the  value  read  from  the  teletypewriter  and  y  is 
the  pressure  in  centimeters  of  water.   Knowing  the  pressure 
readings,  velocities  and  pressure  coefficients  could  then 
be  calculated.   Appendix  A  is  an  outline  of  the  procedure 
used  to  find  the  pressure  coefficients. 

C .   TUFTS 

To  visualize  various  flow  patterns,  tufts  of  yarn  were 
taped  over  the  turret,  nosepiece,  and  base  platform.   Turbu- 
lent boundary  layer  flow  in  the  turbulent  separated  region 
could  be  physically  seen  by  the  violent  motions  of  individual 
tufts.   Likewise,  quiescent  flows  were  demonstrated  by  non- 
moving  tufts.   A  comprehensive  flow  mapping  for  various  flow 
conditions  is  presented  in  the  next  chapter. 
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VI.   EXPERIMENTAL  RESULTS 

A.   TEST  PROCEDURE 

1.  No  Suction 

The  turret,  nosepiece,  base  platform,  and  fairing 
were  set  up  on  the  tunnel  floor  as  depicted  before  in  Figure 
IV-9.   Care  was  taken  to  ensure  all  tufts  were  capable  of 
full  motion  throughout  the  testing.   The  wind  tunnel  was 
activated  and  a  flow  mapping  was  made  with  no  blower  suction. 
The  data  acquisition  system  was  manually  started  and  pressure 
coefficients  of  the  various  surface  ports  were  calculated 
according  to  Appendix  A. 

2.  Suction  Through  a  2.5  Inch  Diameter  Hole 

Two  2.5  inch  diameter  holes  were  drilled  on  both 
sides  of  the  base  platform  at  the  location  from  Chapter  III. 
After  the  wind  tunnel  was  up  to  full  speed,  the  suction  was 
turned  on.   Flow  patterns  again  were  noted  and  mapped  while 
pressure  data  was  recorded.   The  blower  suction  was  varied 
from  full-closed  to  100%  capacity  through  increments  of  10% 
in  order  to  determine  minimum  suction  for  quiescent  flow. 
Total  quiescent  flow  was  never  achieved. 

3 .  Suction  with  the  Aft  Area  Covered 

To  obtain  better  results,  two  large  L-shaped  holes 
were  cut  out  of  the  plywood  on  either  side  of  the  nosepiece. 
The  holes  could  then  be  fitted  with  the  two-piece  adjustable 
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hole  control  or  a  porous  plate  (as  shown  before  in  Figure 
IV-10) .   For  this  part  of  the  test,  the  aft  portions  of  the 
hole  control  covered  the  area  closest  to  the  nosepiece. 
Suction  was  varied  and  data  was  recorded. 

4 .  Suction  with  the  Forward  Area  Covered 

The  same  test  as  in  part  3  above  was  conducted,  but 
the  forward  portion  of  the  hole  was  now  the  covered  area. 

5 .  Suction  Through  a  2  Inch  Slot 

Both  pieces  were  then  screwed  to  the  base  platform. 
A  two  inch  slot  was  left  directly  in  the  center  to  determine 
suction  performance.   This  test  simulated  moving  a  vortex 
"hole"  out  from  the  nosepiece  area  from  where  it  was  posi- 
tioned previously  in  part  1.   The  location  of  the  slot  was 
varied  in  the  general  area  with  very  little  noticeable  flow 
changes. 

6 .  Suction  Through  the  Porous  Plate 

Two  porous  aluminum  plates  were  screwed  over  both 
base  platform  holes  and  similar  tests  were  run.   The  porous 
plates  were  used  to  determine  whether  they  would  help  or 
hinder  the  suction  performance. 

7 .  Suction  Through  Large  Uncovered  Holes 

All  covering  plates  were  removed  during  the  final 
testing.   This  proved  to  be  the  most  optimum  condition  during 
the  experimental  research. 
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B.   FLOW  MAPPING  RESULTS 

To  visualize  the  test  flow  patterns,  figures  were  drawn 
using  the  tufts  as  guides.   Tufts  demonstrated  three  basic 
types  of  motion:   smooth  flow,  mild  turbulence,  and  very- 
strong  vorticity. 

It  was  noticed  during  testing  that  turbulence  would  de- 
crease when  the  suction  was  increased  in  steps.   All  cases 
examined  never  fully  reached  total  quiescent  flow. 

1 .  No  Suction 

Figure  VI- 1  is  a  flow  plot  without  suction  through 
the  base  platform.   Notice  the  extreme  turbulence  depicted 
by  the  vorticity  arrows,  which  is  typical  of  fully  separated 
flow.   There  is  an  upflow  from  the  270  degree  point  at  the 
base  to  the  210  degree  point  close  to  the  top  of  the  turret. 
The  nosepiece  is  also  directing  airflow  in  an  opposite  di- 
rection downwards.   The  entire  nosepiece  and  aft  turret  show 
disturbed  flow  at  almost  all  locations. 

2 .  Suction  Through  a  2.5  Inch  Diameter  Hole 

A  much  quieter  flow  pattern  can  be  seen  along  the 
nosepiece  surface  in  Figure  VI-2.   Notice  that  the  separa- 
tion point  has  moved  further  aft.   However,  the  lower  portion 
of  the  turret  is  still  very  turbulent  with  opposing  flow 
vectors.   The  results  of  this  experiment  forced  the  cutting 
of  a  larger  suction  area  on  both  sides  of  the  base  platform. 
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3 .  Suction  with  the  Aft  Area  Covered 

Figure  VI-3  illustrates  the  improved  performance  with 
a  much  larger  suction  hole  in  the  base  platform.   Notice 
the  aft  hole  area  closest  to  the  nosepiece  is  covered  up  by 
one  part  of  the  two-piece  adjustable  aluminum  hole  control. 
The  turbulent  vorticity  at  the  turret  base  has  ceased,  but 
the  flow  is  now  more  unsteady  nearer  the  top  as  compared  to 
Figure  VI-2.   The  airflow  at  the  bottom  section  of  the  nose- 
piece  is  forced  towards  the  suction  hole,  causing  increased 
disruption.   Flow  along  the  base  platform  is  also  being 
sucked  through  the  hole. 

4 .  Suction  with  the  Forward  Area  Covered 

When  the  forward  portion  of  the  suction  hole  is 
covered  (Figure  VI-4) ,  there  is  a  noticeable  delay  in  separa- 
tion around  the  lower  base  of  the  turret.   Most  of  the  tufts 
are  lying  flat  in  the  direction  of  the  suction  hole.   There 
is  still  some  mild  turbulence  coming  over  the  top  aft  por- 
tion of  the  turret  and  downwards  in  an  "S"  shaped  path. 

5 .  Suction  Through  a  2  Inch  Slot 

As  noted  before,  this  test  was  to  simulate  moving 
the  vortex  hole  out  from  the  nosepiece  body.   From  Figure 
VI-5,  it  can  be  seen  that  the  amount  of  airflow  being  drawn 
off  is  not  enough  to  control  separation.   Various  position- 
ing of  the  flow  control  plates  proved  fruitless  for  quies- 
cent flow  patterns. 
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6 .  Suction  Through  the  Porous  Plate 

Figure  VI-6  depicts  the  porous  plate  covering  the 
entire  suction  area.   Notice  that  the  upper  aft  portions  of 
the  turret  still  show  mild  turbulence  while  the  separation 
point  was  moved  to  about  the  225  degree  position.   The  porous 
plate  had  a  solidity  of  about  60%.   It  is  interesting  to 
note  that  the  larger  hole  areas  perform  better  even  when  the 
flow  is  reduced  by  a  fixture  such  as  a  porous  plate. 

7 .  Suction  Through  Large  Uncovered  Holes 

As  was  expected,  the  optimum  performance  was  demon- 
strated with  the  maximum  possible  area  open  for  suction. 
Figure  VI-7  shows  the  best  flow  field  in  all  the  major  runs. 
Notice,  however,  that  total  quiescent  flow  was  never  fully 
achieved  during  these  tests.   Increasing  the  suction  rate 
would  usually  cause  the  turbulence  to  decrease  at  some 
points,  but  the  design  did  not  produce  adequate  performance. 

C.   PRESSURE  COEFFICIENT  RESULTS 

Pressure  coefficient  values  from  section  BB  (Figure  V-l) 
were  calculated  for  various  runs.   The  runs  included:   no 
suction,  50%  suction  through  2.5  inch  diameter  holes,  10% 
suction  through  large  uncovered  holes,  and  50%  suction 
through  large  uncovered  holes.   Pressure  coefficient  values 
from  section  BB  were  plotted  because  that  section  was  the 
major  region  of  turbulence  for  all  test  runs.   Results  were 
plotted  against  the  theoretical  surface  pressure  distribution 
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for  a  cylinder.   Pressure  taps  at  the  180  degree  position 
were  covered  by  the  nosepiece,  therefore  those  coefficients 
were  not  plotted.   Figure  VI-8  is  a  graph  of  the  pressure 
coefficients  for  the  various  test  runs  compared  to  theory. 
Notice  how  increased  suction  and  increased  suction  hole  size 
demonstrate  a  pressure  distribution  close  to  the  theoretical 
plot. 

Table  VI-1  gives  a  pressure  coefficient  distribution  of 
the  nine  nosepiece  ports  during  test  runs.   Recall  that  this 
area  usually  showed  very  calm  flow  conditions. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  use  of  the  nosepiece  design  with  suction  generally- 
improved  flow  conditions  over  the  laser  turret.   However, 
total  quiescent  flow  was  never  fully  achieved  during  all 
test  cases.   Increased  suction  port  area  in  the  base  plat- 
form permitted  better  flow  performance. 

The  nosepiece  design  was  not  capable  of  stabilizing  shed 
vortices.   Most  test  cases  with  suction  demonstrated  smooth 
flow  over  the  top  of  the  turret. 

Based  on  the  calculations  conducted  during  design,  the 
computer  model  was  not  an  effective  tool  for  flow  prediction 
The  nosepiece  design  requires  more  accurate  information  for 
construction  due  to  its  complexity. 

B.  RECOMMENDATIONS 

The  idea  of  stabilized  vortices  for  blunt  body  flow  con- 
trol should  not  be  discounted.   A  more  adequate  three  di- 
mensional computer  model  is  needed  to  fulfill  design 
requirements.   A  vorticity  measuring  device  would  help  to 
determine  strength  and  location  of  trapped  vortices. 
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APPENDIX  A 
EVALUATION  OF  THE  PRESSURE  COEFFICIENT 

Another  form  of  turret  flow  performance  can  be  determined 
by  using  the  following  equation 


2f-^=  -  f  (A-l) 

Ad 
where  — *■   is  called  the  pressure  coefficient  (C  ) .   The  value 

q  p 

p   is  the  static  pressure  at  the  point  of  interest  measured 
from  the  turret  or  nosepiece  taps.   The  static  pressure  in 
the  wind  tunnel  is  p^  and  the  freestream  dynamic  pressure  is 
q.   The  difference  between  wind  tunnel  total  pressure  (p  ) 
and  wind  tunnel  static  pressure  (p^)  determine  q.   In  other 
words,  equation  (A-l)  can  be  expressed  as 


An     "s     "oo 


Since  the  pressure  coefficient  is  non-dimensional,  values 
obtained  from  the  teletype  can  be  directly  entered  into 
equation  (A-2)'  without  using  a  "centimeters  of  water" 
conversion. 
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TABLE  III-l 

COORDINATES  OF  INITIAL  PANEL  DESIGN 

(From  Figure  III-5) 

POINT  NUMBER  X,Y  COORDINATES 

1  (0,  0) 

2  (0.134,  0.5) 

3  (0.5,  0.866) 

4  (1.0,  1.0) 

5  (1.5,  0.866) 

6  (2.119,  0.625) 

7  (2.387,  0.625) 

8  (2.452,  0.417) 

9  (2.774,  0.286) 

10  (2.940,  0.327) 

11  (3.071,  0.4762) 

12  (3.131,  0.625) 

13  (3.25,  0.625) 

14  (3.488,  0.625) 

15  -   (8.964,  0.625) 

16  (11.375,  0) 

17  (8.964,  -0.625) 
13  (3.488,  -0.625) 
19  (3.25,  -0.625) 
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TABLE  III-l  (Continued) 


POINT  NUMBER  X,Y  COORDINATES 

20  (3.131,  -0.625) 

21  (3.071,  -0.4762) 

22  (2.940,  -0.327) 

23  (2.774,  -0.286) 

24  (2.452,  -0.417) 

25  (2.387,  -0.625) 

26  (2.119,  -0.625) 

27  (1.5,  -0.866) 

28  (1,  -1) 

29  (0.5,  -0.866) 

30  (0.134,  -0.5) 
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TABLE  III-2 

COORDINATES  OF  NOSEPIECE  VARIATIONS 

(From  Figure  III-6) 


POINT  NUMBER 


X,Y  COORDINATES 


6 

7 

8 

9 

10 

11 

12 

20 

21 

22 

23 

24 

25 

26 


(2.119 
(2.446 
(2.773 
(2.982 
(3.101 
(3.191 
(3.190 
(3.190 
(3.191 
(3.101 
(2.982 
(2.773 
(2.446 
(2.119 


0.119) 

0.149) 

0.179) 

0.25) 

0.327) 

0.446) 

0.625) 

-0.625) 

-0.446) 

-0.327) 

-0.25) 

-0.179) 

-0.149) 

-0.119) 


46 


TABLE  III-2  (Continued) 
(From  Figure  III-7) 


POINT  NUMBER 


X,Y  COORDINATES 


6 

7 

8 

9 

10 

11 

12 

20 

21 

22 

23 

24 

25 

26 


(2.119 
(2.625 
(2.893 
(3.042 
(3.131 
(3.160 
(2.982 
(2.982 
(3.160 
(3.131 
(3.042 
(2.893 
(2.625 
(2.119 


0.298) 

0.179) 

0.119) 

0.155) 

0.250) 

0.381) 

0.607) 

-0.607) 

-0.381) 

-0.250) 

-0.155) 

-0.119) 

-0.179) 

-0.298) 


(From  Figure  III-8) 


POINT  NUMBER 


X,Y  COORDINATES 


(2.00,  0.369) 
(2.357,  0.357) 
(2.625,  0.149) 
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TABLE  III-2  (Continued) 
(From  Figure  III-8) 


POINT  NUMBER       .     X,Y  COORDINATES 

9  (2.833,  0.167) 

10  (3.012,  0.286) 

11  (3.131,  0.446) 

12  (3.190,  0.625) 

20  (3.190,  -0.625) 

21  (3.131,  -0.446) 

22  (3.012,  -0.286) 

23  (2.833,  -0.167) 

24  (2.625,  -0.149) 

25  (2.357,  -0.357) 

26  (2.00,  -0.369) 
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TABLE  III-3 
COORDINATES  OF  FINAL  PANEL  DESIGN 

POINT  NUMBER  X,Y  COORDINATES 

1  (0,0) 

2  (0.134,  0.5) 

3  (0.5,  0.866) 

4  (1.0,  1.0) 

5  (1.5,  0.866) 

6  (2.119,  0.298) 

7  (2.30,  0.179) 

8  (2.55,  0.119) 

9  (2.65,  0.155) 

10  (2.73,  0.250) 

11  (2.75,  0.381) 

12  (2.68,  0.607) 

13  (2.83,  0.625) 

14  (3.070,  0.625) 

15  (8.547,  0.625) 

16  (10.958,  0.0) 

17  (8,547,  -0.625) 

18  (3.070,  -0.625) 

19  (2.83,  -0.625) 

20  (2.68,  -0.607) 
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TABLE  III-3  (Continued) 


POINT  NUMBER  X,Y  COORDINATES 

21  (2.75,  -0.381) 

22  (2.73,  -0.250) 

23  (2.65,  -0.155) 

24  (2.55,  -0.119) 

25  (2.30,  -0.179) 

26  (2.119,  -0.298) 

27  (1.5,  -0.866) 

28  (1.0,  -1.0) 

29  (0.5,  -0.866) 

30  (0.134,  -0.5) 
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TABLE  III-4 
TANGENTIAL  VELOCITIES  AND  PRESSURE 
COEFFICIENTS  OF  FINAL  PANEL  DESIGN 


PANEL 
(between  points) 

TANGENTIAL 
VELOCITIES 

PRESSURE 
COEFFICIENTS 

1 

,    2 

0.495 

0.755 

2 

,    3 

1.339 

-0.794 

3 

,  4 

1.783 

-2.178 

4 

,  5 

1.650 

-1.723 

5 

,  6 

0.787 

0.380 

6 

,    1 

-0.230 

0.947 

7 , 

,    8 

-0.544 

0.704 

8, 

,  9 

-0.194 

0.963 

9, 

,  10 

-0.081 

0.994 

10 

,  11 

-0.145 

0.979 

11, 

,  12 

-0.0 

1.00 

12, 

,  13 

1.347 

-0.814 

13. 

,  14 

1.117 

-0.249 

14, 

-  15 

1.092 

-0.192 

15, 

,  16 

1.037 

-0.074 

16, 

,  17 

-1.037 

-0.074 

17, 

18 

-1.092 

-0.192 

18 

,  19 

-1.117 

-0.249 
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TABLE    III-4     (Continued) 


PANEL 

TANGENTIAL 

PRESSURE 

(between 

points) 

VELOCITIES 

COEFFICIENTS 

19, 

20 

-1.347 

-0.814 

20, 

21 

0.0 

1.00 

21, 

22 

-0.145 

0.979 

22, 

23 

-0.081 

0.994 

23, 

24 

-0.194 

0.963 

24, 

25 

-0.544 

0.704 

25, 

26 

-0.230 

0.947 

26, 

27 

-0.787 

0.380 

27, 

28 

-1.650 

-1.723 

28, 

29 

-1.783 

-2.178 

29, 

30 

-1.339 

-0.794 

30, 

1 

-0.495 

0.755 
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TABLE  VI -1 
PRESSURE  COEFFICIENTS  OF  NOSEPIECE 
(Tap  Locations  from  Figure  V-2) 


PRESSURE 

TAP 

NO 

SUCTION 

10 

%  SUCTION 

50%  SUCTION 

LOCATION 

SUCTION 

2.5"  HOLE 

LARGE  HOLE 

LARGE  HOLE 

1 

0.55 

1.68 

0.17 

-0.20 

2 

0.36 

0.62 

0.35 

0.06 

3 

0.58 

0.72 

0.52 

0.29 

4 

-0.34 

-0.63 

-0.67 

-1.67 

5 

-0.13 

0.18 

-0.52 

-1.75 

6 

-0.54 

0.31 

-0.33 

-1.49 

7 

0.37 

-3.04 

-3.21 

-6.71 

8 

0.41 

-2.13 

-2.14 

-4.62 

9 

0.38 

-2.10 

-1.36 

-3.81 
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Streamlines 


Laser  Turret 


Fairing 


Suction  Port 


FIGURE  1-1 


TRAPPED  VORTEX/SUCTION  FAIRING, 
FROM  CRAIG  [Ref.  6] 
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rn^<0 
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(a)   Non  Viscous  Flow 


(b)   Laminar  Flow 


(c)   Turbulent  Flow 


FIGURE  II-l.   FROM  PAO  [Ref.  5],  STREAMLINE  FLOW  AND 

DISTRIBUTION  FOR  A  CYLINDER 
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FIGURE  II-2. 


VELOCITY  DISTRIBUTION  OF 
A  FORCED  VORTEX 
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FIGURE  II-3. 


VELOCITY  DISTRIBUTION  OF 
A  FREE  VORTEX 


Free  Vortex 


Forced  Vortex 


FIGURE  II-4 


VELOCITY  DISTRIBUTION  OF  A 
RANKINE  COMBINED  VORTEX 


56 


t 


5- 


BURGERS-  ROTT 


RANK  IN  E 


FIGURE    II-5. 


NORMALIZED  GRAPH  OF  THE 
BURGERS -ROTT  AND  RANKINE 
VELOCITY  DISTRIBUTIONS 
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FIGURE  III-6.   FIRST  NOSEPIECE  DESIGN  VARIATION 
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FIGURE  III-7.   SECOND  NOSEPIECE  DESIGN  VARIATION 
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FIGURE  III-8.   THIRD  NOSEPIECE  DESIGN  VARIATION 
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FIGURE  IV- 2. 


PLENUM  CHAMBER  AND  FAIRING 
(NOSEPIECE  AND  TURRET  REMOVED 
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FIGURE  IV- 3.   BLOWER  AND  DUCTWORK,  BACK  SIDE 
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FIGURE  IV-4.   FILLER  PIECE  CONSTRUCTION 
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FIGURE  IV-5.   NOSEPIECE,  TOP  VIEW 
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FIGURE  V-l.   TURRET  PRESSURE  TAP  LOCATIONS 
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FIGURE  V-2.   NOSSPIECE  PRESSURE  TAP  LOCATIONS 
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CGMMON/VAL/X(50) tY(50) ,X 1 , Y 1 , THETA( 5G ) , B ( 50) ,£(50  J 

COMMON/ DA/PI, I,X2,Y2,F,PAI 

COMMON/ CA1/NN, I I, ANGLE, ALPHA, R, PHI, C I, C2,PHIJ,PHIJ1 

INTEGER    N,I tJ  ,K,M,N1 

REAL  NNt lit  ANGLE t ALPHA fRf PHI, C li C2fPHIJ» PHI  J  1 

REAL  A( 50,50)  ,0(50,50) 

REAL  ERR,EE, STORE, C(  50) , COLO (50)  ,  V< 50 ) ,CP (50  I 

M=99 
C  ENTER  THE  NUMBER,  N  OF  PANELS  -  MODIFY  200  FORMAT  TO  SUIT 

N=30 

Nl-N+1 

NN=FLOAT(N) 

PI=ATAN( 1.01/45.0 

PAI=PI*180. 

PAIM=-PAI 

READd,*)  ALPHA, ERR, ICASE 
C  ALPHA  SPECIFIES  THE  LCCATICN  OF  START  OF  FIRST  PANEL 
C  ERR  SPECIFIES  THE  REQUIRED  ACCURACY  FOR  ITERATION 

REA0(1,*)  X2,Y2,F 

DO  1  1=1,  Nl 

II=FLCAT(I) 
C    COORDINATES    OF    START    OF    PANEL    I    TO    6E    FED    IN 

IFdCASE.EQ.DCALL    DATAi 

IF(ICASE.EQ.2)CALL    DATA2 

1  CONTINUE 
DO    2    1=1, N 

IFdCASE.EQ.DCALL  CASE1 
IF(ICASE.EQ.2)CALL  CASE2 
WRITE(6,210)     X(I) ,Y( I) ,THETA(I) 

210    FQRMAT(1H    ,3(iX,F6.3)) 

2  CONTINUE 
DC    3    1=1, N 

Xl  =  (X(I)+X(H-lJ)/2.0 
Yl  =  (Yd  )+Y(I+l))/2.0 
DO    4    J=1,N 

R=UX1-X(  J+l) )**2+<Yl-Y( J+l) )**2)/( (Xl-XU) )  **2 
1         +(Y1-Y(J) )**2) 
PHIJ=ATAN2(Y1-Y< J),X1-X(J)) 
PHIJ1=ATAN2(Y1-Y(J+1)  ,X1-X( J+l) ) 
C  IF(Y( J) .GT.Y1.AND.Y( J+l) .Le. Yl ) PHIJ= 8 .0*ATAN ( 1.0) +PHIJ 

PHI=PHIJ-PHIJ1 

IF(PHI.GT.PAI  )PHI=PHI-2.*PAI 
IF(PHI.LT.PAIM)PHI  =  PHH-2.*PAI 
ANGLE=THETA(J)-THETA(I) 

A(  I, J)=-PHI*CCS( ANGLE) -0.5*ALOG(R)*SIN( ANGLE) 
IF(I.EQ.J)     A(  I,J)=4.0*ATAN(1.0) 
D(  I,J)  =  PHI*SIMANGLE)-0.5*AL0G(R  J*COS(ANGLE) 
4    CONTINUE 

3  CONTINUE 

C    ENSURE    THE    NUMBER    AFTER    THE    ,    IS    THE    NUMBER    OF    PANELS 

200    F0RMAT(1H    ,  3 (  10F6 .3/ ) , / ) 
C    WE    NOW    WISH    TO    SOLVE    FCR    THE    PANEL    SOURCE    STRENGTHS    C(i) 
C    MODIFY    NEXT    SECTICN    FOR    A    BETTER    METHOD    CF    GETTING    C(I) 
00    10    I=1,N 
C(I)=C.C 
10    CONTINUE 

DO    20    K=1,M 
EE=0.0 
DO    30    1=1, N 
COLDd)=C(I) 
STCRE  =  B(  I) 
DO    40    J=1,N 

IF( J.NE.I)STOPE=STCRE-A(I,J)*C( J) 
40    CONTINUE 

C(  I)=STCRE/A( 1,1) 
EE=EE+A8S(C(I )-CCLD< I) ) 
30    CONTINUE 
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IFIEE.LT.ERR)GG    TC    50 
20    CONTINUE 
50    CONTINUE 

HRITE(6,220)K 
220    F0RMATC1H    ,23HNUMBER    OF    ITERATIONS    -     ,IZ) 

*RITE(6,250) 
250    FCRMAT<1H    ,29HTHE    PANEL    STRENGTHS,    C(I)     ARE) 
WRITE(6,200) (C(  I) ,1=1, N) 
C    THIS    CONCLUCES    THE    SOLUTION    CF    LINEAR    EQUATIONS    FOR    C(IJ 
C    NEXT    OBTAIN    THE    VELOCITIES    TANGENTIAL    TO    THE    PANELS,    V(I) 
00    60    1=1, N 
V<  I)  =  E(  I) 
DO    70    J=1,N 
V(  I)=V(  I  )+0<  I  ,J)*CU) 
70    CONTINUE 
C    ALSO    OBTAIN    THE    CORRESPONDING    PRESSURE    COEFFICIENTS    CP(I) 

cp(i>=i.o-(v(  n**2j 

60  CONTINUE 

hRITE(6,230) 
230  FORMAT! 1H  ,35HTHE  TANGENTIAL  VELOCITIES,  V(I)  ARE) 

WRITE(6,200)( V( I),I=1,N) 

URITE<6,240) 
240  F0RMAT(1H  ,36HTHE  PRESSURE  COEFFICIENTS,  CP(I)  ARE) 

WRITE(6,200)(CP<  I)  ,1=1, N) 

WRITE(6,*)X2, Y2,F 

STOP 

END 

SUBROUTINE  DATA1 

CCMMGN/VAL/X(50> , Y ( 5  0 ) , XI, Yl , THE TA{ 5C  ) , 6 ( 50)  , E ( 50 ) 

C0MM0N/DA/PI,I,X2,Y2,F,PAI 

COMMON/ DA1/NN, 1 1, ANGLE, ALPHA, R, PHI, C1,C2,PHIJ,  PHI  J  1 

REAL  NN, II,ANGLE,ALPHA,R,PHI,C1,C2,PFIJ,PHIJ1 

ANGLE =(  180.0-ALPHA-3  60.C*UI-1.0)/NNI*PI 

X(I)=COS(ANGLE) 

Y(I)=SIN(ANGL£) 

RETURN 

END 

SUBROUTINE  CASE1 

COMMON/VAL/X( 50) ,Y(50) , XI, Yl, THE TA( 5C  ), B ( 50)  , EI50) 

COMMON/ DA/PI, I ,X2,Y2,F,PAI 

THETAd  )=ATAN2(Y(  1  +  1  )-Y  (  I )  ,  X(  1  +  1  )-X  (  I  ) ) 

B( I)=SIN(THETA(I )  ) 

E( I)=CGS(THETA(I) ) 

RETURN 

END 

SUBROUTINE    DATA2 

COMMON/ VAL/X( 50) ,Y(5C)  , XI, Yl , THE TA( 50  )  ,B < 50)  , E(50) 

COMMON/ 0 A/PI,  I,X2,Y2,F,PAI 

READ(lt*J    X(I),Y(I) 

RETURN 
ENO 

SUBROUTINE  CASE2 

COMMON/ VAL/X( 50) ,Y(50) , XI, Yl , THETA( 50 ) ,B ( 50) ,E(50) 

COMMON/ DA/PI, I,X2,Y2,F,PAI 

REAL  NX, NY 

X1  =  (X(I)+XU+1)  )/2. 

Y1=(Y(I )+Y( I+l) )/2. 

P12=(X1-X2)**2+(Y1-Y2)**2 

R22=< Xl-X2)**2+( Y1+Y2J**2 

U=F/(2*PAI)*< (Y1-Y2)/R12-(Y1+Y2) /R22) 

V=-F/(2*PAI)*((X1-X2)/R12-(X1-X2)/R22) 

THETA(I)=ATAN2(Y(  I+l)-Y<  I) ,X(I  +  1)-X(  I)) 

NX=CCS(THETA(  D-PAI/2.  I 

NY  =  SIN(THETA<  I  )-P A  1/2. J 
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C  B<  I)  =  SIN(THETMI  )) 

B(  I)=U*NX+V*NY+SIN(THETA(I)  ) 

TX=COS(THETA( I)) 

TY=SIN(ThETA< I)) 
C  E<  I)=COS(THETA(U) 

E(  I)  =  U*TX  +  V*TY*CCSUHETA(I)  i 

RETURN 

END 
C 

C       INPUT    DATA    FOLLCfcS 
C 

22.5,0.0001,2 
2.4220,0.3200,0.6 
0.0, 0.0 
0.1340,0.5 
0.5,0.8660 
1.0,1.0 
1.5,0.8660 
2.119,0.2976 
2.3000,0.1786 
2.5500,0.1190 
2.6500,0.1548 
2.7300,0.2500 
2.7500,0.3810 
2.6300,0.6071 
2.83,0.625 
3.0700,0.625 
8.5470,0.625 
10.958,0.0 
8.5470,-0.625 
3.0700,-0.625 
2.83,-0.625 
2.6800,-0.6071 
2.7500,-0.3810 
2.7300,-0.2500 
2.6500,-0.1548 
2.5500,-0.1190 
2.3000,-0.1786 
2.119,-0.2976 
1. 5,-0.8660 
1.0,-1.0 
0.5,-0.8660 
0.1340,-0.5 
0.0,0.0 
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A  silent  16  ram  movie  is  available  on  loan  to  interested 
parties.   The  movie  shows  turret  tuft  motion  with  and  with- 
out flow  control  employed  for  different  nosepiece  designs. 
Requests  should  be  directed  to  Professor  Allen  E.  Fuhs ,  Code 
6  7Fu,  Department  of  Aeronautics,  United  States  Naval  Post- 
graduate School,  Monterey,  California   93940.   Telephone: 
Commerical,  (408)  646-2948,  or  AUTOVON,  878-2948. 
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